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Abstract 

We investigate the deconfining phase transition in QCD4 and QED4 at finite temperature us- 
ing a perturbative deformation of topological quantum field theory (TQFT). A modified maximal 
abelian gauge (MAG) is utilized in the analysis. In this case, we can derive the linear potential 
studying the 2D theory through Parisi-Sourlas (PS) dimensional reduction. The mechanism 
of deconfining phase transition is proposed. It is geometrical to discuss the thermal effect on 
the linear potential. All we have to do is to investigate the behavior of topological objects as 
such instantons and vortices on a cylinder. This is the great advantage of our scenario. This 
mechanism is also applied in the case of QED4. The phase structure at the high temperature 
of QED is investigated using the Coulomb potential on a cylinder. It coincides with the result 
in the lattice compact U (1) gauge theory. Also, QCD with MAG has the property called the 
abelian dominance, which enables us to discuss the deconfinement of QCD4. 
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1 Introduction 



Quark confinement is one of main problems in QCD. This phenomenon is realized at least in 
the low energy region (IR region) or at low temperature. However, due to asymptotic freedom 
the coupling constant becomes large in the IR region, where perturbation theory would not be 
reliable and applicable. Therefore, quark confinement should be explained from non-perturbative 
aspects of QCD. This phenomenon could be explained well using the dual super-conductor 
vacuum scenario based on monopole condensation Q. However, this scenario is not sufficient 
and many ideas have been still proposed. One of the recent scenarios, which explains the quark 
confinement in "continuum" QCD, is based on using a perturbative deformation of topological 
quantum field theory (TQFT) g g,§. This scenario describes the quark confinement very well. 
In particular, if we choose a modified maximal abelian gauge (MAG), which is a kind of partial 
gauge fixings ||, then the linear potential between a static quark and antiquark appears [Q]. 
This means the quark confinement in the Wilson criterion [r|. 

The properties of QCD medium at finite temperature have been the subject of the intense 
study. It undergoes a drastic change as the temperature increases. It is believed that confined 
quarks and gluons are liberated from the certain temperature and the system is deconfined. 
That is, deconfining phase transition should be caused. For high temperature a characteristic 
energy of quarks and gluons traveling through the medium is high, and the effective coupling is 
small. The system presents the quark-gluon plasma (QGP). Since the effective coupling is small, 
we can use the perturbation theory and the perturbative calculation of its self-energy leads to 
the thermal mass of gluons and quarks. These mean the Debye screening effect as in the usual 
plasma. 

In this paper, we study the deconfining phase transition of the finite temperature QCD4 
(i.e., thermal QCD) with MAG. In a previous paper Q, we have investigated the difference of 
the thermal phase transition mechanism between Lorentz type gauge fixing and MAG. In both 
cases the TQFT sector can have the phase structure. In Lorentz type gauge fixing this structure 
is essential to the deconfinement of full QCD4, while in the MAG the phase structure of TQFT 
sector cannot have relevance in full QCD4. Therefore, it has been unclear how we can explain 
the deconfinement phase transition. In this paper, we propose the deconfinement scenario in 
thermal QCD4 with MAG, in which all we have to do is to investigate the behavior of topological 
objects such as instantons and vortices on a cylinder. For example, in an SU(2) QCD4 we need 
to consider the instantons of 2D 0(3) non-linear sigma model (NLSM2 or CP 1 model) on a 
cylinder ||. However, this theory has asymptotic freedom and their instanton solution has the 
size parameter. Therefore the treatment seems rather difficult because the usual dilute gas sum 
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ansatz is reliable in the very restricted region. Hence we concretely argue the deconfmement 
scenario in compact QED4 in this paper. In this case, the TQFT sector becomes an 0(2) NLSM2 
and only to consider their vortices on a cylinder, which have no issues as the above. We also 
can investigate the behavior at high temperature using the Coulomb potential on a cylinder. 
Moreover the result can be applied to the abelian projected effective theory, which is an abelian 
gauge theory with asymptotic freedom. Finally some prospects and issues in studying QCD4 
are also discussed. 

Our paper is organized as follows. In section 2, we review the method of a perturbative 
deformation of the TQFT at zero temperature. We choose a modified MAG, which have some 
additional terms and the gauge fixing term has an OSp{A\2) symmetry. Quark confinement could 
be explained analyzing TQFT sector, in which the non-perturbative information of confinement 
has been encoded. This sector is equivalent to a certain 2D theory through Parisi-Sourlas 
dimensional reduction mechanism ||, which is caused due to an OSp(A\2) symmetry. We mainly 
investigate the case of G = SU(2), in which the TQFT sector becomes an 0(3) NLSM on a 
plane (zero temperature). It has also instanton solutions and the linear potential between 
quark-antiquark pair could be induced by the instanton effect. In section 3 we consider the 
finite temperature case. The evaluation of the rectangular Wilson loop at zero temperature 
become the evaluation of two Polyakov loops' correlator. Therefore, we can investigate the 
phase transition from a viewpoint of 2D instantons on a cylinder. It is expected that the main 
role in the deconfmement phase transition is determined by the behavior of these. In section 4, 
we investigate the phase structure of compact QED4. The method of a perturbative deformation 
can also be applied to QED4 and it is known that zero temperature compact QED4 has confining 
phase at large coupling |K|. In this case, the TQFT sector is an 0(2) NLSM2 on a cylinder. 
Of course, at zero temperature it is an 0(2) NLSM2 on 2-plane, and a confining potential is 
induced by the Coulomb gas of vortices. The phase transition is described by the celebrated 



Berezinskii-Kosterlitz-Thouless (BKT) phase transition |11]. It is natural extension to consider 
the Coulomb gas on a cylinder. Thus we can investigate the high temperature region of QED4. 
We conclude in this case that thermal effect shifts the value of the coupling in which confining 
phase transition is caused. Moreover, we apply the result of compact QED4 to the scenario in 
| |12| by using the abelian dominance. We could explain the deconfmement which is induced by 
the thermal effect. Section 5 is devoted to conclusions and discussions. 
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2 QCD 4 as Perturbative Deformation of TQFT 



2.1 Setup 

In this section, we review the method of a perturbative deformation of TQFT, which is discussed 
m Refs.gg. Firstly, we consider an SU(N) QCD at zero temperature in the 3+1 dimensional 
Minkowskian space-time. The modification in the case of the finite temperature system is 
denoted latter. We don't include matter fields, that is, consider the gluodynamics. The action 
for the gauge group G = SU (N) is 

1 



S = ~2^J d xTx G F ^ r ( 2J ) 
where F^ u is a field strength 

Ffu, = d^A u — d v A^ — i\Ap, A v ] 

and SU (N) generators T A , (A = 1, • • • , iV 2 — 1), which are hermite and traceless, are normalized 

as 

Tt g {T a T b ) = \5 ab . 

In order to construct the quantum field theory from the classical action (|2.1| ), it is adequate 
to use the BRST quantization. Incorporating the (anti-) FP ghost field C{C) and the auxiliary 
field B, we can construct the BRST transformation 5b 

«5 B A M = D„[A]C, 5 B C = iC 2 , 

5 B C = iB, 5b B = 0, (2.2) 

where .D^L4] is the covariant derivative given by 

Dp[A]=d ll -i[A l » ]. 

The gauge fixing term can be constructed from the BRST transformation 5b as 

•Sgf+fp = -i&a J d 4 xG GF+FP [^, C, C, B], (2.3) 

where Ggf+fp[^4^ 5 C, C, B] is determined by the gauge fixing condition . If we perform the 
gauge fixing in the Lorentz gauge in usual manners, then we have 



G GF+FP = Tr G [C(d^ + -B)\ = C A [d^ + -B A j , (2.4) 



where a is a gauge parameter. 
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MAG type gauge fixing 



Lorentz type gauge fixing 



TQFTsector 
2D 0(3) NLSM 



Perturbative 
deformation parts 



U{1) backgrounds 



TQFTsector 
2D 0(4) NLSM 



Perturbative 
deformation parts 



Figure 1: The comparison between Lorentz type gauge and MAG type gauge in G = SU(2). The phase structure 
of the TQFT sector is preserved in the case of the Feynman type gauge fixing. In MAG, however, it is screened 
by the perturbative fluctuation and its information is encoded in the (7(1) backgrounds. 



In this paper, we choose the modified MAGQ 

Ggf+fp^, °i °, B) = ~5 B Tr G/H [A^ + 2iCC] 



-A° 1 .A afl + iC a C a 
2 M 



(2.5) 



where H is the maximal abelian subgroup, which for example is 17(1) 1 in G = SU(N), the 
subscript a denotes non-diagonal generators and 5 B is the anti-BRST transformation: 



& B A IX = D„[A]C 1 5 B C = iB, 
5 B C = iC 2 , 5 B B = 0, B + B = {C,C}. 



(2.6) 



This gauge is MAG with some additional terms. In this case, the gauge fixing term has a 
special symmetry, OSp{4\2) symmetry, which is very powerful to investigate the linear potential 
between a quark-antiquark pair. 

Note that monopoles also exist in MAG. This existence is ensured by the fact that the 
homotopy group ^{G/H) = tt\{H) is non-trivial. These can not exist in the Lorentz type 
gauge. 

Also, modified Lorentz gauge 



G GF+FP LV C, C, B] = 5 b Tt g [ApA" + 2iCC] = 5 B 



l -A A A A » + iC A C A 



is discussed in Ref.||]. The comparison between MAG type and Lorentz type in G = SU(2) is 
shown in Fig.|I[ 

2.2 Decomposition into Topological Trivial and Non-Trivial Sectors 

In this section, we decompose the action of QCD into a topological trivial sector and a non- 
trivial sector which is described by a topological quantum field theory (TQFT), called a TQFT 

*This choice is rather specific, and more general gauge fixing given by 

G GF+FP = Tr G /H [A^ - iaCC] 
is also possible. In this case, the interesting phenomenon, which is discussed in Ref. J^], is caused. 
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sector. This method was proposed in Ref.|| in the context of Lorentz type gauge and extended 
to the MAG in Ref.||. Thus, we may recapitulate QCD as a perturbative deformation from a 
TQFT sector. In this method, the information on the non-perturbative phenomena of QCD, in 
particular quark confinement and phase structure is encoded in the TQFT sector. 
We begin with the following decomposition 

A^ = UV^ + iUd^U\ (2.7) 

where we define as 

% = iUdprf. 

We assume that is given by a finite rotation of Vn- Here J)^ is composed of compact 
degrees of freedom U alone, but UV^U^ is not compact. In below, we assume that the non- 
compact gauge field variable does not have topologically non-trivial configuration and all 
topologically non-trivial configurations come from the compact gauge group variable U alone. 
As suggested in Ref . [p| , it is expected that the perturbation-theoretical study has nothing to 
do with confinement and the real deep reason for the confinement is encoded into the topological 
structure of the gauge group. 

Secondly, we introduce the FP determinant AppLA] defined by 

AfpIA]- 1 J [dU]l[6 (2.8) 

x,A 

where [dU] is the gauge invariant Haar measure and so this determinant is invariant under the 
gauge transformation, 

A FP L4] =A FP [A U ~ 1 ]. 

Then we rewrite the unit as follows, 

i = A FP [A] A^n^(^r) 

x,A 

= AppLO J [dU]\{5 (d^ 1 ) 

x,A 

= A FP [V] I [dU]\{5 {d»Vf) 

J x,A 

= J[dU][d 7 ][d^][d(3]exp 
Here, we define the new BRST transformation <5p as 

£b7 = iP, SbP = 0. (2.10) 



(2.9) 
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By the use of this 5b, eq. (|2.9| ) can be rewritten as 



1 = / [dU][di\[dry][d0\exp 



i \ d x -i5 B G G F+Fp[V^,7,7,/3] 



where Ggf+fpIX^, 7, 7, j3\ is written as 

G G f+fp[^,7,7,/3] = 2Tr G (75%). 
When we insert eq.( |2.1l| ) into the partition function, it is rewritten as 



(2.11) 



(2.12) 



Z[J] = J [dU] [dC] [dC] [dB] J [dV] [d 7 ] [dj] exp 



d 4 * Tr G {F^[V]F^[V}) 
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i5 B G GF+FP [V^ 1 ^,(3}-i6 B G GF+FP [Q^ + UV^U\C,C,B}^ +iSj 



where Sj is a source term given by 



Sj = Id xTr G 



(2.13) 



J» [top + UV^ J + J C C + J c (7 + J B B . 
Note that the transformation law of U and under 5 B and <5b is expressed as 

6 B U = iCU, 5 B U = iCU, 

&BVp = 0, ~5 B V^ = 0. (2.14) 
By the use of these transformation laws, we can rewrite ^bGgf+fp^a 1 + UV^II^] as 
d 4 x<5 B G GF+F p[^ + UV, L U\ C, C, B] 
= J ^bGgf+fpI^, C, C,B] + J d*x (v A M A ^[U] + \v A V B "K AB [U]\ , (2.15) 
where M A [U], and /C AB [J7] are defined as 

M A [U] = S B 5 B [[UT A U ] ) a 9, a ^j , JC AB [U] = 5 B 5 B [{UT A U^) a {UT B U ] ) a ^ 
Finally, we obtain the following expression of the partition function 



Z[J] 



where Stqft is defined by 



[dU] [dC] [dC] [dB] exp (zStqft + iW[U; J M ] 
+ i J d A x(J^ + J C C + J C C + J b b) 



Stqft = — ^bCgf+fp[^, C, C, B] 

I d A x5 B 5 B Ti G/H [fi^fi" + 2iCC] 



(2.16) 



(2.17) 



Here, W[U ; J M ] in eq. (|2.16|) denotes a perturbative deformation and is denned as 
exp (iW[U; J"]) = J [dV] [dj] [<Zy] [dp] exp (iS pQC n [^,7, 7, 0] 

+i J d*x (v^ - \v A V B »K. AB [U]\ \ , (2.18) 
where is the new source term redefined as 

jf = rtj£u-Mi[u], 

and /C describes the interaction between pQCD and TQFT sectors. pQCD denotes the pertur- 
bative QCD (topological trivial sector). Here, the action £ p qcd is defined by 

S p qcd[V^7,j,P] = J d 4 x f-^F^F^iV] - iS B G GF+FP [V^l,P] \ ■ (2.19) 

Note that A4„[i7] and 1C AB [U] are interactions between pQCD and TQFT sector. 

The perturbative deformation W[i7; J^"} is the generating functional of the connected Green 
function of Va, which describes the perturbative deformation sector. This should be calculated 
by the use of the ordinary perturbation theory with the expansion of the coupling constant g 



W[U; J M ] = In ^exp (i J d A x[V A J A >* - V A V B »K 



i 

,2 



ABi 

pQCD 



= 9 Y J d 4 xj d% <WK(y))pQ CD {J^)J v {y) - 5\x - y)^K AB [U\) 

+ higher order of g. (2.20) 

Therefore, W[£/; J^] is expressed as a power series in the coupling constant g and goes to zero 
as g — > 0. It turns out that the full QCD is reduced to the TQFT sector in the vanishing limit 
of coupling constant. Thus we can interpret the term W[£/; J^ 1 ] as the deformation from the 
TQFT sector. 

2.3 Relation between Expectation Values 

Let us discuss below in the Euclidean metric. We can define the expectation value in each sector 
as 

(d . . . O m ) TQFT = J[dU][dC][dC}0 1 ...O m e W {-S TQFT [U,C,C,B}), (2.21) 

(0i...On)pQGD = / [dV}[d^[dj}[d(3}O l ...O n exp(-S pQ cB[V^^,0]), (2.22) 

and reconstruct the expectation value of the full QCD4. If the inserted operator is decomposed 
as f(A) = g(V, U)h(U), then the full expectation value of this is expressed as 

if (A)) QCD — ^)) p qcd ^(^)) TQFT — ((sO 7 ) ^^^^tqft) pQCD • (2.23) 
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T (Imaginary time direction) 



R 



Contour C 



Figure 2: The rectangular Wilson loop. 



Thus we can obtain the full expectation value through the above expectation values in each 
sector. Which of expression eq. fl2.23l ) should be utilized in calculating the expectation values 
depends on the case. In fact, the decomposition of the inserted operator is rather difficult in 
non-abelian gauge group. Our purpose is to investigate the quark confinement and so would like 
to evaluate a non-abelian Wilson loop 

W c = TrPexp (ie j> dx^ 



where the contour C is the rectangular loop as shown in Fig.||. In this case, the expectation 
value of non-abelian Wilson loop is mathematically decomposed as the above by the use of 
non-abelian Stokes theorem [14, 15]. When we consider the case of G = SU{2) for simplicity^, 
it could be expressed as 



(W C [A\) 



QCD 




ie f dx^Tr [asUVf.W 

c 



exp 



pQCD 



3 <t> dx^l 



(2.24) 



TQFT 



where = 2Tr(T 3 M ) and d/i is the invariant Haar measure of the coset space SU(2)/U(1). 
When we expand eq. (|2.24|) perturbatively, it becomes 

„2 



(Wc[A]) QCD = I d/i(C7)^exp 



ie 



dx^m 



c 



TQFT 



r dx** j> dy u D^{x-y) x 



x / dji(U) ^exp 



+ higher order of e. 



ie & dx^W 



c 



2Tr{T 3 U(x)T A U\x)}2Tr{T 3 U(y)T A W (y)} 



TQFT 

(2.25) 



where we used the following relations 



pQCD 



6 AB 6. 



I IV 



5 AB D^(x-y). 



pQCD 4n 2 \x-y\ 



'For G — SU(N), the expression is rather complicated. 
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Here e is a charge of the external source and proportional to g, for example e = g/2 for the 
fundamental representation. Hence, the above expansion is about the power of g. As we will 
see later, the evaluation of the "abelian" Wilson loop 



exp 



(2.26) 



leads to the linear potential. Note that eq.( [2.25| ) does not imply the abelian dominance, which 
is well known feature in the MAG, but tell us what we should evaluate. 

However, in the case of QED4 the decomposition as the above is completely done as follows 

If 



(W C [A]) 

QED — 

(WcM) 

TQFT 

(Wc[V}) pV{1) (2.27) 

where 

A tl = v fl + -ud„u\ n„ = -ud^. 
9 9 

This is similar with the decomposition of the partition function 

Z = Zi nst ■ -ZpU(i) 

in the result of Polyakov's work [jH]], in which the linear potential is derived from Z- mst though 
the theory is on 3D Euclidean. 

2.4 PS Dimensional Reduction to 2D Theory 

Parisi-Sourlas mechanism can dimensionally reduce 4D TQFT sector to 2D theory |J. This is 
because we have chosen the special gauge fixing that has an 05p(4|2) symmetry. Though the 
2D space can be taken arbitrarily, we should take a 2-plane at zero temperature and a cylinder 
at finite temperature as 2D space respectively in order to evaluate Wilson loop to derive the 
linear potential and investigate quark confinement . As a result, the action of TQFT sector 
becomes as follows, 

2ir f 

Stqft = J d 2 xTr G/H [^Q" + 2iCC] . (2.28) 

It describes a coset G/H chiral model on 2D space. 

In the case of SU(2), this action can be rewritten as the action of 0(3) NLSM, 

Stqft = 4 / d2xd v n ■ d M n, n • n = 1. (2.29) 
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Here we used the Euler angle representation of an SU (2) matrix 

U(x) 

exp(|(v? + x))cos| exp(-|(^-x))sin| 



exp ( ixix)^) exp (iO(x) exp (V(z)y 



- exp - X)) sin f exp (-§ 
and parameterized an unit length vector field as 



n(x) 



(<P + X)) 



cos ■ 





n l (x) ^ 




( 


sin^(x) cos ip(x) 


\ 




n 2 (x) 






sin#(x) sin</?(cc) 




V 


n 3 (x) ) 




\ 


cos 9{x) 


) 



(2.30) 



(2.31) 



We can investigate a TQFT sector through an NLSM2. In particular, the expectation values 
of both theories are given as follows, 



{0\ . . . Cn)TQFT 4 — (^1 • • • ^")nLSM 2 ■ 



(2.32) 



2.5 Confinement and Static Potential 

Here, we can evaluate the abelian Wilson loop through 0(3) NLSM instantons as follows, 



{W c [n]) TQFTi = ^exp 



yexp 



2-7ri ( - ) Qnlsm 2 



= ( exp 
TQFT 4 \ 



tqft 4 

2-7ri ( - ) Qnlsm 2 
\9 



(2.33) 



nlsm 2 



where e is a charge of an external source, and 



<3nlsm 2 



d 2 xe^ u n(x) ■ (9 M n(x) x d u n(x)), 



(2.34) 



is an instanton density of an NLSM2. Thus, we can calculate this abelian Wilson loop by the 
use of dilute gas approximation^, and the result is given as 



(Wcp}) TQFTi = exp {-a A) 



(2.35) 



Here A = RT is the area spanned by the contour C, and a, which is a string tension of confining 
string, is given by 



a = 2Be 



-Si 



8tt 

&1 — "Ti 



(2.36) 



*We could improve the instanton calculation including the interaction between (anti-)instantons fl7|| . 
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for the fundamental representation. In eq.( |2.36| ) Si is the 1-instanton action and B is the 
constant derived from the integration of the instanton moduli. When the contribution of the 
perturbative deformation part is included, the full Wilson loop expectation value is written as 



(Wc[A]) QC D = e 



tRT 



3\ e 2 



(2.37) 



and the full static potential between a pair of quark and antiquark is expressed by 



V(R) = aR-(^^f(R) + ..., (2.38) 

where f(R) is a certain function that behaves as f(R) — > 1, (R — > 0). 

Thus the linear potential is induced by the instanton effect in an 0(3) NLSM2 in the leading. 
This means quark confinement in the Wilson criterion. Hence we must consider the behavior 
of the instantons in order to study the confinement and deconfinement essentially. Note that 
the NLSM2 instantons should be interpreted as the points that monopoles' current lines pierce 
the 2D space which has been chosen in dimensionally reducing TQFT4 sector (for details, see 
Ref.Q), and so we can say that monopoles are essentially relevant to the confinement. 

3 Confining Phase and Deconfining Phase Transition 

3.1 Finite Temperature 

Now we would like to consider the finite temperature system coupled to the thermal bath. 



The imaginary time formalism and real time formalism |18|, IS, ^(J are well known procedure 
to investigate a thermal field dynamics (TFD). In both cases, gauge fields obey the boundary 
conditions 

A M (-i/3,x) = A M (0,x) (3.1) 

for an imaginary time direction. 

The twisted boundary conditions for the gauge group element U(x) 

Bi : U{-iP,x) =f/(0,x)e 2 ^/ JV , (/ = 0, ■ • • , N - 1) (3.2) 

can be imposed using the element of the center. For I = it is a periodic one. 
The FP determinant at finite temperature is modified as follows, 

1 7V_1 f 

1 = A ^]^E / [dU]HS^A u ^\x)). (3.3) 
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Due to the thermal effect it is decomposed into N independent sectors. In each sector the gauge 
transformation obeys each boundary condition given by eq. (|3.2| ). It is likely to consider that 
such decomposition should have something to do with the domain wall, which is related to the 
spontaneous discrete symmetry break down, such as the center symmetry Zjy pil l- This point 
remains to be unclear. 

At finite temperature, the expectation value is decomposed to the sum as 

</(^)>qcd - ^ E (i9(V, U)) pQCD h(U))' =nH ( W W^)> tQft) • (3-4) 
3.2 Boundary Conditions in Reduced 2D Theory 

Let us consider the case of N = 2 concretely. In N = 2 using PS dimensional reduction both 
TQFT4 "* (i = 0, 1) sectors are described by the field n(x) obeying the periodic condition for 
imaginary time direction M. 

Boundary conditions of n(x) are derived from the following useful relations 

n A {x)T A = U\x)T 3 U(x), n A (x) = 2Tr G [U {x)T A U ] {x)T% {A = 1, 2, 3). (3.5) 

We find that 11(3;) is invariant under £7(1) transformation generated by T 3 and it can be rotated 



by generators associated with the coset SU(2)/U(1). Also, eq.( 2.2S| ) has the following global 
5C/(2)l symmetry| 

U — >Uh, V/i£ SU{2) L . (3.6) 

Then n(x) transforms as 

n A T A _^ n A(tf T A h j = n' A T A (3.7) 
and we easily see the action of h on n A 



n A — > n' A 



ad(/i t ) A B n B . (3.8) 

B=l 

This means that n(x) should be transformed under 50(3) rotation but it is invariant by an 
action of the center Z2 of SU(2). By the use of eq,(|3.8;), boundary conditions for U(x) can be 
translated into that on field n(x) 

3 

n A (-i(3,a) = ad (5 f ) V B (0,^), ff£Z 2 

B=l 

= n A (0,a) (3.9) 
''In the Lorentz type gauge fixing, TQFT sector becomes 2D 0(4) NLSM2, which has a global chiral symmetry 
S[/(2)l ® SU(2)-r. While in MAG SU(2)l symmetry only exists. Also, this system does not have instanton 
solutions. 
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where a is a spatial coordinate of 2D space. Thus, n(x) obeys a periodic condition. 

In conclusion, boundary conditions of U(x) is irrelevant in MAG from the viewpoint of 
reduced 2D theory and its contribution to the full expectation value is same in each sector, i.e., 
, i 

</(^)>qcd = 2£(fcW W0>§Jft) qcd — > ({g(V,U)h(U)) TQFT ) (3.10) 

i=0 

This point is different from the case of the Lorentz type where the field variables on a reduced 
2D theory obey twisted boundary conditions in each TQFT sector. 

In the case of N > 3, there are some possibilities to take the coset space, so it is nontrivial 
whether similar result is concluded. At least if we take the flag space Fn = SU(N)/U(1) N ~ 1 
(maximal abelian) or the complex projective space C-P^ -1 = SU(N)/(SU(N — 1) x U(l)) as 
the coset, then the similar result seems to be followed [|l5|]. Due to the U(l) factor in the coset 



the structure such as eq.(3.10) would be followed. 



Comment on phase structure of TQFT sector If we consider only TQFT sector without 
perturbative deformation part (or consider pure gauge configuration only), the twist factor 
becomes arbitrary element of SU(2) instead of the center, and more general boundary conditions 
of U(x) are allowed as follows, 

B g : U(-i0,x) = U(O,x)g, (V g € G). (3.11) 

Also, n(x) can obey twisted boundary conditions. In this case, TQFT sector can have the phase 
structure though spontaneously symmetry breaking (SSB), though it is forbidden in 2D theory 
by the novel Coleman- Mermin- Wagner theorem |2^] . This is shown in Ref . |J by calculating the 
effective potential. 

3.3 Polyakov Loop and Confinement at Finite Temperature 

Here, using the imaginary time formalism we consider an SU{N) gauge theory at finite temper- 
ature, in which the order parameters for deconfining phase transition are the expectation values 
of the Wilson lines wrapping k times around the compact time dimensions r = — ix°, 

•k/3 



Pk (x) = TrP exp 



' / oItA t (t, x) 
Jo 



(3.12) 



where P denotes the path-ordering and P)-(x), (k 6 Z) are also called Polyakov loops. The 
theory has a "Ln symmetry, which allows us to impose the twisted boundary conditions 

U(0,x) = e 27rli/N U(0,x) (Z = (),-■■ ,JV-1) (3.13) 
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for gauge group elements. Under this gauge transformation the Polyakov loop transforms as 

P k (g)^e 2 * m ' N P k {x). (3.14) 

Therefore (P k (x)) = means that Z^r symmetry is unbroken. A Polyakov loop corresponds to 
an external quark source (in the fundamental representation of SU(N)). The free energy of the 
system (heat bath) is increased by adding such a source. If we write this additional free energy 
as F q , the expectation value of a Polyakov loop can be expressed as 

(P k (x))~e-P F «. (3.15) 

From eq.( |3.15| ), {P k (x}) = implies that the free energy cost is infinite and an isolated quark 
can not exist in the theory. While if (P k (x)) ^ 0, it is finite and an isolated quark has finite 
energy, i.e., the theory no longer confines. In summary, 

(P k (x)) = ==> confining phase and Z^r symmetry is unbroken, 
(P k (x)) ^ ==> deconfining phase and Z^r symmetry is broken. 

Thus Zjv symmetry has been considered to be related to deconfinement transition. Also, Z^r 
symmetry is important to decide the order of the phase transition. 

In general, QCD4 is believed to be confined at low temperature. Therefore, the correlator of 
two Polyakov loops is expected to behave as 

(P fe (P)P_ fe (0)) ~ 0, (R = \R\ — ► 00), (3.16) 

F q q = aR, (a : string tension) 



at low temperature. That is, it should show the exponentially decay law. Also, eq,( 3.16p implies 
|(P)| 2 = 0, and so confinement. 

Therefore we would like to evaluate the correlator of the Polyakov loops in order to study the 
confinement. Note that the rectangular Wilson loop at zero temperature becomes the Polyakov 
loops' correlator at finite temperature as shown in Fig.||. Thus, 

(Pl(£)P_i(0)> = {WcM) (3.17) 

is followed. Therefore, in the same way as the case of the zero temperature, we can derive the 
linear potential and obtain the following result 

<Pi(i?)P_i(0)) ~e^ CTi? . (3.18) 



14 



R 




Pi 



A 



7^ 



Area A(C) = iZT 



This part is canceled. 
Correlator of Polyakov loops 



Figure 3: The correlator of Polyakov loops. The expectation value of a Wilson loop W is equivalent to a 
correlator of the Polyakov loops Pi and P-\. 

Note that string tension a depends on the temperature T as 

a = 2B(T)e~ s \ B(T) ~ 1/T, S 1 = (3.19) 

because the instanton moduli integral is restricted on the cylinder and the effective coupling 
depends on the temperature. Therefore the string tension depends on the temperature contin- 
uously, and decrease as the temperature increases. 

Also, a perturbative deformation part would be modified by the thermal effect. In particular, 
we can derive the Yukawa type potential, which means the Debye screening effects as in the usual 
plasma instead of the Coulomb potential at the sufficiently high temperature region. In this time, 
we use the standard calculation in TFD and use the hard thermal loop (HTL) approximation 



[23]. Therefore, calculating the contribution of the gluon self-energy to its propagator at 1-loop 



level with HTL approximation, we obtains the 

x xab 1 

<W^G/)) P QCD = j-P e~ m ^-y\, m D = -g 2 T 2 C A , (3.20) 

p 4tt\x — y\ o 

where mo is the thermal gluon mass and Ca = N is the group constant. 
3.4 Deconfinement Phase Transition 

The result eq.( |3.18| ) is favorable at least in low temperature region but at high temperature 
region the system should be deconfined. Therefore it is expected that the string tension should 
vanish above a certain temperature, i.e., deconfining phase transition should be caused by the 
thermal effect^]. This problem has been remarked in our previous paper 0. Let us recall that the 
■ Unlike in the Lorentz type gauge fixing, the phase structure of TQFT sector would not be retained in MAG 
once perturbative deformation part is included. Therefore the deconfining phase transition must be explained in 
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linear potential is induced by topological objects in our scenario. In order for the linear potential 
to vanish, the effect of topological objects has to be ignored. For example, a pair of topological 
objects needs to form the bound state and behaves as "neutral molecules", or topological objects 
decay. In such case, the linear potential would vanish since it is induced by non-zero topological 
charge in the rectangular Wilson loop. From the viewpoint of reduced 2D theory, the thermal 
effect is realized through the radius of the cylinder. While at low temperature the cylinder is 
expanded infinitely and behaves as 2-plane, at high temperature cylinder collapses and behaves 
as ID line at sufficiently high temperature. This implies that it is so simple and geometric to 
study the thermal effect on the linear potential and the deconfining phase transition. All we 
have to do is to investigate the behavior of topological objects on a cylinder. This is the great 
advantage of our deconfinement scenario. 

Thermal effect on the instantons was discussed in Ref . p4|| . At high temperature the cylinder 
radius constrains the instanton size and the behavior of these is changed as the temperature 
increases. Large-size instantons are suppressed and so instanton effects can be reliably calculated 
at sufficiently high temperature. That is, the dilute gas approximation is rather valid, and so 
the interaction of instantons should be suppressed. In fact, at sufficiently high temperature the 
instanton effect could be ignored. 

However, it is difficult to deal with the instanton beyond the dilute gas approximation and 
to calculate concretely. Therefore we consider the case of a compact U(l) gauge theory i.e., 
compact QED without matter. In QCD, we must deal with instanton and so worry with the 
various problem like infrared problem. But there is no problem in QED as such. If the gauge 
group U(l) is compact, abelian monopoles and the confining phase exist. In this case, TQFT 
sector is an 0{2) NLSM on the cylinder and the linear potential is induced by the vortices of 
an 0{2) NLSM, which should be interpreted as abelian monopoles from 4 dimensional view. 
Vortices are realized by the globally neutral Coulomb gas on the cylinder (on 2-plane at zero 
temperature). It was showed in Ref. J10[| that the confining-deconfining phase transition at zero 
temperature can be described by BKT phase transition at certain critical coupling. Using the 
expression of the propagator on the cylinder we could find the thermal effect on this phase 
transition. In the next section, we will discuss in detail. 

Order of Phase Transition Here, we comment on the order of the deconfining phase transi- 
tion. There are two possibilities of the deconfinement mechanism. One is that the string tension 
vanishes discontinuously because of the vanishing of topological object effect as is discussed 
different way from in the case of the Lorentz type. 
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Figure 4: Thermal effect is realized as circumference of a cylinder or width of a strip. Conformal transformation 
w = ^ lnz maps the z-plane onto the strip. 



above. This corresponds to the first order phase transition. Thus, the both possibilities of the 
deconfining phase transition seems to exist. Another is that the string tension continuously 
vanishes at sufficiently high temperature because of e~ Sl dumping or the T-dependence of B. 
This case would correspond to the second order transition. 

In general, it is considered that the order of deconfinement phase transition depends on the 
gauge group. The deconfinement of 4D SU(N) QCD without quarks can be related to the 
ferromagnetic-paramagnetic transition of "Lm spin model with the ferromagnetic interactions 
and Zjv symmetry plays a crucial role in the confinement problem. In particular, in the case of 
N = 2, it is related to 3D Z2 spin model (Ising), and the deconfinement transition is the second 
order. Also, in the case of N = 3 |2f|, it is related to 3D three-state Potts model, the phase 
transition of which is the first order one also expects that a deconfinement phase transition in 
an SU(3) theory, is of first order. 

In the above scenario, it is unclear how we can explain the order of deconfinement transition. 
We could possibly explain it from the intense study of the behavior of instantons on reduced 2D 
theory gl. 



4 Deconfining Phase Transition in Compact QED 4 



The phase structure of compact QED4 at zero temperature has been worked out in Ref.|lC]. It 
has confining phase above the certain value of the coupling. Whether such a phase exists or not 
is decided due to the compactness of U(l). If U(l) is non-compact, such phase can not exist. 
In next subsection, we study the thermal effect on the confining phase. 
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Phase structure of 2D Coulomb gas 

Critical temperature 
1/8tt 

I 1 



Dipole phase Plasma phase with Debye screening 

Dipolcs dissociate. 



Figure 5: 2D Coulomb gas has two different phases. Over the critical temperature l/8n, the system is the 
plasma phase with Debye screening, and therefore mass gap exists. Below 1/87T it is the dipoles phase, in which 
Coulomb charges form dipoles. The system has a long range correlation and no mass gap. 



4-1 Compact QED4 at Finite Temperature 

We consider the case of compact QED4, which has the confining phase in the strong coupling 
region (UV region). In the similar way, its TQFT4 sector becomes an 0(2) NLSM2, vortices of 
which induce the confining potential. The contribution of vortices to the partition function is 
described by 2D classical Coulomb gas. One of the difference between QCD and QED is that 
there exists no scale parameter (like A parameter in QCD) in QED and also no size parameter 
in vortex solutions of 0(2) NLSM2. Hence the treatment is rather simple and strict. When we 
consider the finite temperature case, the TQFT sector becomes 0(2) NLSM2 on the cylinder 
and vortices contribution could be described by the Coulomb gas on the cylinder. 

The phase structure of 2D Coulomb gas is well known as shown Fig.||. Here, we would like 
to consider the Coulomb gas on the cylinder. The propagator on the cylinder is given by 



G(w 



w 



Am 

47T 



e 



2 + 2^ Re(u7 + ^ 



(4.1) 



where w = a + ir and r is the imaginary time. We can find eq.( |4.1[ ) by constructing this propa- 
gator on the complex plane and then mapping it on a cylinder using conformal transformation 
0. Then, requiring that G(w — w') is a function of w — w', an additional term appears. The 
propagator eq.(|4.1|) is rewritten as 



G(w - w') 



Air 



In 



2 1 

+ — Re(ui — w') 



In the low temperature region ((3 — > 00), 



G(w 



1 

47T 
1 



In 



2tt 
T 



(w 



2/5 



(4.2) 



In \w 



w 



(4.3) 



Another derivation of (Kit) is given in Appendix A. 
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As is expected, this eq.(|4.2j) becomes usual 2D Coulomb potential, and behaves as the propagator 
on the plane. The constant term can be ignored because of the neutrality of the Coulomb gas. 
Next, we consider in the high temperature region (f3 — > 0). In the case of Re(w; — w') > 0, 



G{w — w') 

In the case of Re(w — w') < 0, 

G(w - w') 



In 



1 

47T 



w 



(4.4) 



-±-ln\-l\ 2 + ±Re(w- W >), 



1 



Re(w — w'). 



Therefore we find in the high temperature region ((3 — > 0) 

1 



G(w - w') 



■ t \Re(w — w') 



1 



a — a 



(4.5) 



(4.6) 



2/5 1 n 2/3' 

This is ID Coulomb potential. As is expected, the propagator behaves as on a line. Note 
that the factor 1/(3 appears. This factor does not appear if we naively start from ID Coulomb 
gas. The origin consists in the finite cylinder radius. This factor is very important in the later 
discussion. 

Let us discuss the behavior of Coulomb gas. It behaves as ID Coulomb gas and its partition 
function is given by 



'nc 



oo 

^ (n!)< 

n=0 v ; 



J d(Tj exp 



i>3 



C 



-S(i) 



(4.7) 



where is a single vortex action. Note that the temperature of this Coulomb gas system is 
9 = [3g 2 = g 2 /T. The ID Coulomb gas is exactly solvable [27]. Its behavior in large 6 and 
small 9 is well known. Its phase structure is very similar with the 2D Coulomb gas. For small 
9 the ID Coulomb gas behaves as the gas of free "molecules", made up of H — charges pairs 
bound together. On the contrary, for large 9, the charges are completely deconfined, forming an 
electrically neutral "plasma" of 2n free particles. Now we are considering the high temperature 
region {(3 — > 0, that is T — > oo) i.e., the small 9, and the ID Coulomb gas behaves as a gas of free 
"molecules". This implies that the linear potential always vanishes irrelevantly to the definite 
coupling constant if the temperature is sufficiently large. Of course, if the coupling g 2 becomes 
much larger than T, that is we consider such an energy scale where g 2 becomes too large, the 
theory is confined again. The relative measurement between the coupling g 2 and the physical 
temperature T is important to decide whether confining or deconfining. 
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Confining 



Figure 6: Phase diagram of the lattice compact U(l) gauge theory in 4D. The horizontal line is the coupling 
constant and the vertical one is the temperature. 

As a result, we could say that compact QED4 is deconfined at the sufficiently high temper- 



ature. This nicely corresponds to the lattice compact U(l) gauge theory |25|, |28|]. Its phase 
diagram is shown in Fig.^. Therefore we conclude that our scenario describes the confinement- 
deconfinement in compact QED4 at zero and finite temperature very well. 

4-2 Abelian Dominance and Deconfining Transition 

In this section, using the previous results and the abelian dominance, we discuss the deconfining 
phase transition of an SU(N) QCD4. 

In MAG, it is expected that off-diagonal gluons are massive and diagonal gluons are massless 
and that the diagonal component dominate in a sense of Wilsonian effective action. It is con- 



firmed by Monte Carlo simulation on a lattice that off-diagonal gluons are massive 30]. As 
an analytical derivation has been given at least in the TQFT sector based on the PS dimensional 
reduction to the coset G/H NLSM2 0. Once this fact would be assumed, we could deal an 
SU(N) QCD4 as an U (l)^ -1 abelian gauge theory in the low energy region, which is N — 1 
copies of compact QED4. That is, the low energy effective theory becomes an abelian gauge the- 
ory. Therefore we can apply the result of QED4 in previous subsection. This scenario has been 
proposed in Ref . p^j , which also can be applied to the finite temperature case. This 
gauge theory which is derived under the assumption of abelian dominance, has asymptotic free- 
dom in contrast to usual (compact) QED4. Hence the coupling constant becomes large in the IR 
region as the same as QCD4. It is clear that the thermal effect could cause the deconfinement 
phase transition at sufficiently high temperature from the result of previous subsection if we be- 
gin with this low energy effective theory. Hence in this scenario our deconfinement mechanism 
works well. 
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Figure 7: 2D Coulomb gas is intimately related to 2D sine-Gordon model and 2D massive Thirring model. 



In the case of N = 3 the critical temperature has been estimated in Ref.[36], based on the 
dual Ginzburg-Landau description. 

5 Conclusions and Discussions 

We have investigated deconfining transition of QCD4 and QED4 at finite temperature. In our 
quark confinement scenario, using PS reduction the analysis is reduced to 2D theory on the 
cylinder, in which the thermal effect is realized through the radius of cylinder. Therefore, it 
is very simple and geometrical to investigate thermal effect on the deconfining transition. In 
order to discuss the deconfinement it is enough to study the behavior of topological objects on 
the cylinder with several radii. In particular, we could investigate the phase structure of the 
compact U(l) gauge theory in " continuum" . This result agree with the result of lattice. 

Taking high temperature limit is equivalent to compactification on a small circle S 1 , and so 
we could discuss the relation between different dimensions, for example QCD4 and QCD5. 5D 
theory has been very interesting since it is widely discussed in the context of theory with extra 
dimensions. It is interesting work to investigate the connection between 5D and 4D QCD with 
the procedure in this paper. Also, our method might be applied to the system with Higgs fields 
such as Georgi-Glashow model. In Ref.[31], in which the Georgi-Glashow model in 3D at finite 



temperature is discussed, at high temperature the system become 2D theory and BKT phase 
transition appears. 

It is interesting to consider more general case G = SU(N) rather than G = SU(2) that 
is mainly discussed in this paper. If N > 3, there is the possibility of choosing the coset. 
The simplest example is a maximal torus SU(N)/U(l) N ~ 1 . In this case, N — 1 kinds of 
monopoles exist. However, a single monopole is enough to induce the linear potential, and 
so SU(N)/(SU(N - 1) x C/(l)) ^ CP N ~ l is also possible. The confinement needs at least 
one U(l) factor in the above scenario. In CP N_1 case, the similar argument with the case of 
G = SU(2) seems to be possible [[19]. The large N behavior of CP^ -1 model is well under- 
stood. Therefore we might investigate the large N behavior of QCD4 with our method of the 
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perturbative deformation. 

In this paper, we investigated the phase structure of the compact U(l) gauge theory at 
sufficiently high temperature and strong coupling using the behavior of the ID Coulomb gas. 
The intimate relation among the classical Coulomb gas, massive Thirring model and sine-Gordon 
model at zero temperature is well known as indicated in Fig.f?]. This relation at finite temperature 



is discussed in Refs.[|32;, |33(| . In particular, we could discuss the phase structure of the compact 



U(l) gauge theory in terms of sine-Gordon model on the cylinder instead of studying the classical 



Coulomb gas behavior [34]. The phase structure of sine-Gordon model at finite temperature is 
discussed in the Ref . []35|] . Also, it is attractive to include dynamical fermions and consider a 
finite chemical potential, chiral symmetry or flavor quantum number iVf. 

The relation to the result of the dual Ginzburg-Landau theory is also interesting. In the case 
of N = 3 the critical temperature of the deconfinement has been estimated in Ref. [j36[l , based on 
the dual Ginzburg-Landau description. 
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Appendix 

A The propagator of a periodic boson on a cylinder 

The action 

i rP 



S = - J dx I dxi d^dfd (A.l) 



describes massless free boson on a cylinder. 
The propagator is given by 



n=—oo 
-oo 



i ^ r dh e - iHx -y) r dse -s(k^) ) (A 2) 

P n=-oo •* ^ ^° 
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where k 2 = + k\, ko = 2im/[} (n € Z) is Matsubara frequency. In the above calculation we 
introduced the mass term [i in order to avoid the infrared divergence. The above expression 



( A. 2 ), using the modified Bessel function Ko(z), is rewritten as 

D(x -y) = — ^ K o (W (xo -Vo- n(i) 2 + (x[ - y{) 2 

n=—oo 

In the limit fj, — > 0, D(x — y) becomes 

- +00 

D(x-y) ~ -— ^ In (jiy/ (xq - y - n(i) 2 + (x[ - yi) 2 ) , 
1 



(A.3) 



^-ln ( pt/3Wcosh f^O^i -2/1 



( ^"( x o - 2/o) 



Here, we introduce the complex coordinates 

w = xi + ix , {w = xi- ix ), w' = yi + iy , (w' = yi - iy ), 



(A.4) 



(A.5) 



and rewrite D(x — y) by the use of these coordinates (AJ5). Simple calculation leads to the 
following expression, 



D(x-y) 



47T 



In 



e 1 3 — e P 



2 1 
+ — Re(w + w') 



1 



(A.6) 



2/3 " v " ' " ' 1- 

The last term in the eq.( |A.6| ) implies the infrared divergence. But this term can be ignored 
because of the neutral condition of Coulomb gas. Thus, we obtained the propagator ([4.1|), Also, 
we can relate the sine-Gordon model on the cylinder with the classical Coulomb gas on the 
cylinder, using the propagator ( [4 . 1[ ) . 
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